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The dynamical creation of entanglement between three-level atoms coupled to the common vacuum is inves- 
tigated. For the class of bound entangled initial states we show that the dynamics of closely separated atoms 
generates stationary distillable entanglement of asymptotic states. We also find that the effect of delayed sud- 
den birth of distillable entanglement occurs in the case of atoms separated by a distance comparable with the 
radiation wavelength. 
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I. INTRODUCTION 

Dynamical creation of entanglement by the indirect inter- 
action between otherwise decoupled systems has been studied 
by many researchers mainly in the case of two-level atoms 
interacting with the common vacuum. When the two atoms 
are separated by a distance small compared to the radiation 
wave length A,, there is a substantial probability that a pho- 
ton emitted by one atom will be absorbed by the other and 
the resulting process of photon exchange produces correla- 
tions between the atoms. Such correlations may cause that ini- 
tially separable states become entangled (see e.g. lUl-H]). The 
case of three-level atoms is very interesting for many reasons. 
First of all, in a system of coupled multi - level atoms hav- 
ing closely lying energy states and interacting with the vac- 
uum, quantum interference between different radiative transi- 
tions can occur, resulting in coherences in a system which are 
known as vacuum - induced coherences. Beside the usual ef- 
fects such as collective damping and dipole-dipole interaction 
involving non-orthogonal transition dipole moments Slalj ra- 
diative coupling can produce here a new interference effect in 
the spontaneous emission. This effect manifests by the cross 
coupling between radiative transitions with orthogonal dipole 
moments |7] and is strongly dependent on the relative orien- 
tation of the atoms |8, 9]. All such collective properties of 
the system influence the quantum dynamics, which can sig- 
nificantly differ from a corresponding single atom dynamics. 
On the other hand, the theory of entanglement between the 
pairs of such atoms is much more complex then in the case of 
qubits. As is well known, there is no simple necessary and suf- 
ficient condition of entanglement, since the Peres-Horodecki 
separability criterion Ill0l[ll]l only shows that the states which 
are not positive after partial transposition (NPPT states) are 
entangled. But there can exist entangled states which are pos- 
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itive after this operation and such states are not distillable yj 
Hence all entangled states can be divided into two classes, one 
contains free entangled states that can be distilled using lo- 
cal operations and classical communication (LOCC) and the 
other consist of bound entangled states for which no LOCC 
strategy can be used to extract pure state entanglement. Since 
many effects in quantum information exploit pure maximally 
entangled states, only distillable states can be directly used 
for quantum communication. Contrary, nondistillable states 
states involve some kind of irreversibility: we need pure en- 
tanglement to create them, but no pure entanglement can be 
obtained back from them 1 13]. 

The dynamics of entanglement between three-level atoms 
with vacuum-induced coherences was previously investigated 
mainly in the context of creation or degradation of NPPT 
states ll 14 [l5ll . In particular it was shown that for small dis- 
tance between the atoms the system decays to a stationary 
state which can be entangled, even if the initial state was sep- 
arable [140. On the other hand, if the distance is compara- 
ble to the radiation wavelength, the dynamics brings all initial 
states into the asymptotic state in which both atoms are in 
their ground states but still there can be some transient entan- 
glement between the atoms lUSll . 

In the present paper, we investigate the dynamics of bound 
entangled initial states. There are some results concerning 
decoherence and disentanglement of bound entangled states 
ilql . but here we focus on the process of dynamical creation 
of distillable entanglement due to the collective damping and 
cross coupling between the three-level atoms. For the spe- 
cific bound entangled initial state and small interatomic dis- 
tance we show that the asymptotic state is both entangled and 
distillable. Thus we obtain a stationary free entanglement. 
The same result is also valid for other initial states includ- 
ing separable states. So this dynamics of three - level atoms 
distinguishes distillable states, since all nontrivial asymptotic 
entangled states are also distillable. For larger distances, the 
dynamics of the bound entangled initial state is very peculiar: 
the system very quickly disentangle and only after some finite 



time there suddenly appears a distillable entanglement. (The 
similar phenomenon of delayed sudden buth of entanglement 
was observed in the case of two-level atoms UM ). So also 
in this situation the physical process of spontaneous emission 
can create some transient distillable entanglement out of the 
initially prepared bound entanglement. 



II. MIXED - STATE ENTANGLEMENT AND 
DISTILLATION 

A. Distillability of entanglement 

Distillability of mixed entangled state p is the property that 
enables to convert n copies of p into less number of k copies 
of maximally entangled pure state by means of LOCC [18'\. 
It is known that all pure entangled states can be reversibly 
distilled |19] and any mixed two-qubit entangled state is also 
distillable 1,20,1 . In general case, the following necessary and 
sufficient condition for entanglement distillation was shown 
in Ref. III2II : the state p is distillable if and only if there exists 
n such that p is n - copy distillable i.e. p®" can be filtered to 
a two-qubit entangled state. This condition is however hard to 
apply, since conclusions based on a few copies may be mis- 
leading 1I21I1 . More practical but not necessary condition is 
based on the reduction criterion of separability 112211 . The cri- 
terion can be stated as follows: if a bipartite state p of a com- 
pound system AB is separable, then 



where 



PA(8)l-p>0 and l(g)pB-p>0 



PA = trBP, PB=trAP 



(II. 1) 



As was shown in Ref. ll23ll . any state that violates dll.ll ) is dis- 
tillable, so if 



PA(8)ll-p^0 or ]l(g)pB-p^O 



(11.2) 



the state p can be distilled. The condition ill.lj is easy to 
check and we will use it in our discussion of dynamical as- 
pects of distillability. 



B. Peres-Horodecki criterion and bound entanglement 

To detect entangled states of two qutrits, we apply Peres- 
Horodecki criterion of separability |10, 11]. From this crite- 
rion follows that any state p for which its partial transposition 
p^^ is non - positive (NPPT state), is entangled. One defines 
also negativity of the state p as 



A^(p) 



Ip'^lltr- 



1 



(11.3) 



A^(p) is equal to the absolute value of the sum of the nega- 
tive eigenvalues of p^^ and is an entanglement monotone ll24ll . 
however it cannot detect entangled states which are positive 



under partial transposition (PPT states). Such states exist f25|] 
and as was shown in Ref. 112J, are not distillable. They are 
called bound entangled PPT states. Up to now, it is not known 
if there exist bound entangled NPPT states [26]. 

To detect some of bound entangled PPT states we can use 
the realignment criterion of separability ll27[ 12811 . The crite- 
rion states that for any separable state p of a compound sys- 
tem, the matrix R{p) with elements 



I (g) {fi\R{p)\n) (g) |v) = {m\ (g) {n\ p \^i) (g) |v) 



(11.4) 



has a trace norm not greater then 1 . So if the realignment 
negativity defined by 



Nr{p) =max 0, 



l^(p)l|tr-l 



(11.5) 



is greater then zero, the state p is entangled. In the case of two 
qubits, the measure dll.Sb cannot detect all NPPT states 112911 . 
but for larger dimension the criterion is capable of detecting 
some bound entangled PPT states ll27ll . 



III. TIME EVOLUTION OF THREE-LEVEL ATOMS 

To study the dynamics of entanglement between three-level 
atoms we consider the model introduced by Agarwal and Pat- 
naik |7]. We start with the short description of the model. 
Consider two identical three - level atoms ( A and B) in the V 
configuration. The atoms have two near - degenerate excited 
states |la), |2a) (a = A,B) and ground states |3a). Assume 
that the atoms interact with the common vacuum and that tran- 
sition dipole moments of atom A are parallel to the transition 
dipole moments of atom B. Due to this interaction, the pro- 
cess of spontaneous emission from two excited levels to the 
ground state take place in each individual atom but a direct 
transition between excited levels is not possible. Moreover, 
the coupling between two atoms can be produced by the ex- 
change of the photons, but in such atomic system there is also 
possible the radiative process in which atom A in the excited 
state I 1a) loses its excitation which in turn excites atom B to 
the state |2b). This effect manifests by the cross coupling be- 
tween radiation transitions with orthogonal dipole moments. 
The evolution of this atomic system can be described by the 
following master equation |7] 



dp 

dt 



i[H,p] + {L'^+L'^ 



-L^')P 



(III.l) 



where 



H^Y^^kiiot^oi + afioi, 



3k) 



k=l 



+ £ n,,(a?3a3i« + a?2Cfi'3«) 



(III.2) 



a=A,B 



and for a =A,B 



L«p = Y^ y,3 (2af,po£ - a'^.af.p po^o?,) (III.3) 



k=l 



moreover 



was found in Ref. lll4ll . It is given by 



L^^p 



^ ^ rB(2a?,pa^3«-a^3X^p-paeXJ 

k=la=A,B 

^ r,.,(2ofipo23" - a23"o?iP - pa23Xi 



(111.4) 



-A, o"^ is the transition operator from |^a) to \Ja) 



In the equations (IIIOI i. jlini l and (|nL4] l. 

and B for a = 

and the coefficient j^ represents the single atom spontaneous 

- decay rate from the state \j) ( j = 1,2 ) to the state |3). 

Since the states |la) and |2a) are closely lying, the transition 

frequencies CO13 and CO23 satisfy 

CO13 ~ CO23 = COq 
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na is A for a = B where 



1 



(III.5) 



J^ = ^ (P22 + 2p33 + P44 + 2p77 - 2Re P24 - 4Rep37) 

Z = T (P36 - P38 - P76 + P78) 

>V = - (P26 + P28 + 2p39 - P46 - P48 - 2p79) 

y=Z (P22 + P44 + 2p66 + 2p88 - 2Rep24 - 4Rep68) 



Similarly, the spontaneous - decay rates 



Yi3 ~ Y23 = Y 

The coefficients Y j-^ and Q.j^ are related to the coupling be- 
tween two atoms and are the collective damping and the dipole 
- dipole interaction potential, respectively. The coherence 
terms Yy^ and Hyc are cross coupling coefficients, which cou- 
ple a pair of orthogonal dipoles. Detailed analysis shows the 
cross coupling between two atoms strongly depend on the rel- 
ative orientation of the atoms and there are such configura- 
tions of the atomic system that Y^c — ^vc — and the other 
configurations for which Y^c 7^ 0, i2,,f 7^ 0. Moreover, all the 
coupling coefficients are small for large distance R between 
the atoms and tend to zero for /? — > 00. On the other hand, 
when 7? — > 0, i2i3, ii23 and H,,^ diverge, whereas 

ri3,r23^Y, and r„,^o 



The time evolution of the initial state of the system is 
given by the semi-group {r,},>o of completely positive map- 
pings acting on density matrices, generated by the hamiltonian 
(IIII.2b and the dissipative part L^ +LP + L^^. The properties 
of the semi-group crucially depends on the distance R between 
the atoms and the geometry of the system. Irrespective to the 
geometry, when R is large compared to the radiation wave- 
length, the semi-group {r,},>o is uniquely relaxing with the 
asymptotic state |3a) ® |3b). Thus, for any initial state its en- 
tanglement asymptotically approaches 0. But still there can be 
some transient entanglement between the atoms. On the other 
hand, in the strong correlation regime (when R — > 0), the semi- 
group is not uniquely relaxing and the asymptotic stationary 
states are non-trivial and depend on initial conditions. The ex- 
plicit form of the asymptotic state Pas for any initial state p 
with matrix elements p^./ (with respect to the canonical basis) 



V = T (-P23 - P27 + P43 + P47 + 2p69 - 2p89) 



(III.6) 



and 



t=\-2x-2y 



Depending on the initial state, the asymptotic state ( IIII.5I I can 
be separable or entangled. In the next section we will study 
distillability of Pas for some initial states. 



IV. GENERATION OF STATIONARY DISTILLABLE 
ENTANGLEMENT 

As was shown in |14], the negativity of the asymptotic 
states ( |III.5t can be obtained analytically in the case of di- 
agonal (i.e. separable) initial states. For such states, only the 
parameters jc, y and t are non-zero and the asymptotic negativ- 
ity reads 



A^(pas) = 2 



\/4(^'+/)- 



(IV. 1) 



Note that every nontrivial asymptotic states from that class is 
entangled. Now we show that this entanglement is free i.e. all 
such asymptotic states are distillable. To do this, we show that 
Pas corresponding to diagonal initial states, violate reduction 
criterion ( III. 11 1. Indeed, since for such states 

/.=c 0\ 
0x0000000 
000000x00 
OOO3/OOOOO 

tTfiPas «) 1 - Pas = OOOOjOOOO (IV.2) 

OOOOOOOyO 

OOxOOOaOO 

OOOOOjO/70 

VO 0000000 c/ 



where 



l-2x- 



-y, 



-2y, 



and the matrix on the right hand side of ( IIV.2l i has two negative 
leading principal minors (other minors are positive), so 



Similarly 



trBPas«)l-Pas ^0 



1 (g) tlA Pas - pas ^ 



The interesting examples of nontrivial asymptotic states are 
given by the separable initial states where the one atom is in 
the excited state and the other is in the ground state or two 
atoms are in different excited states. In all such cases, the 
created entanglement is free and can be distilled. 

Now we consider the possibility of creating free stationary 
entanglement from the bound initial entanglement. As the ini- 
tial states we take the family II30I1 



Pa = y|n)(^0| + fP+ + ^f- 



3<a<4 (IV.3) 



where 






(^1 



1a>®|2b> ■ 



and 



(^|2a>®|1b>- 



'^|2A>®|3fi>+^|3^>®|lB>j 



-^|3a>®|2b>+^|1a>®|3b>; 



The states ( IIV.3I ) are constructed as follows: we prepare the 
maximally entangled state |^o) and add some specific noise 
resulting in the mixing of |M^o) with separable states f+ and 
P . For a special choice of mixing parameter such prepared 
states have positive partial transposition but are entangled, as 
can be shown by computing the realignment negativity. For 
Pa it is given by 



A^R(Pa) 



1 
21 



V3a2- 15a +19-1 



(IV.4) 



and is obviously positive for 3 < a < 4. So the states ( |IV.3t 
are bound entangled and the entanglement initially present in 
l^o) cannot be extracted from them, for any number of copies 
of the states. It is worth to notice that recently the bound en- 
tanglement was created experimentally in the system of three 
qubits [31]. 

Although the states ( IIV.3l l are not diagonal, one can check 
that the corresponding asymptotic states have the same form 
as in the diagonal case. In fact, for all initial states pa there is 
only one asymptotic state pas given by 



x — y— — and t — — 
^ 56 14 



By the above discussion, this state is entangled and moreover 
its entanglement is distillable. Thus we have shown that the 
physical process of spontaneous emission in the radiatively 
coupled three - level atoms can transform initial bound entan- 
glement into free distillable entanglement of the asymptotic 
state. 



V. DELAYED CREATION OF DISTILLABLE 
ENTANGLEMENT 

In this section we study in details the evolution of entangle- 
ment of the bound entangled initial states (IIV.3l l for 3 < a < 4, 
beyond the strong correlation regime. In that case, the asymp- 
totic state is trivial, but some transient entanglement between 
the atoms can be produced. For simplicity, we consider such 
atomic configuration for which the cross coupling coefficients 
are equal to zero. One can check that the initial states ( IIV.3I ) 
will evolve into the states of the form 



/Pll 











Pl5 











Pl9\ 





P22 





























P33 











P37 

















P44 

















P51 











P55 











P59 

















P66 





P68 











P73 











P77 























P86 





P88 





\P91 











P95 











P99/ 



Pa(0 



(V.l) 
where all non-zero matrix elements are time dependent. 

Numerical analysis indicates that during the time evolution 
the realignment negativity ill-Si of the initial state very rapidly 
goes to zero, so the system almost immediately disentangle. 
To consider possible creation of free entanglement, let us first 
check if PPT condition can be violated during such evolution. 
After taking the partial transposition, the state dV.ll ) becomes 



Pa(t 



,PT 



/pii 

P22 Pi5 
P33 
P51 p44 






VP37 



P37\ 


pi9 


P55 P68 

P66 P59 

P91 P77 

P95 P88 

P86 P99/ 

(V.2) 
One can check that determinant d of the matrix ( |V.2| l equals 








d = (P22P44 - |pl5|^)(p33p77 - |P19P)(P66P88 - |P59p) 



X (pilP55P99-P55|P 



37 



-pll|P86|^) 



(V.3) 



We can show numerically that ( IV.3I ) changes the sign, since 
the last factor is positive for all t <tN and becomes negative if 
t > In, for some t^ > 0, and the remaining factors are positive. 
Moreover, all other leading principal minors of the matrix 
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FIG. 1 : The time evolution of the factor F = p 1 1 psg pgg — P55 1 P37 1 — 
PlllPsel^ in l lV.3t for the initial state ( IIV.3b with a = 3.6 and the 
interatomic distance R = 0.2 A.. 




FIG. 2: The time evolution of the factors G = rT,^r-n — IP37I (dotted 
line) and H = rs^rgg — Ipssl (solid line) in l lV.6t for the initial state 
l [rv3] > with a = 3.6 and R = 0.2A,. 



(IV.2I 1 are always positive. So the evolution of the bound entan- 
gled state (IIV.3I 1 has the interesting property: for all t <tN the 
states Pa(f) are PPT and then suddenly they become NPPT 
states (see FIG. 1). 

Now we discuss distillability of the states Pa(f)- Since 
we cannot exclude the possibility that there are NPPT states 
which are non-distillable, we try to apply the reduction cri- 
terion of entanglement. As we know from the discussion in 
Sect. II, any state violating this criterion is necessarily distill- 
able. By direct computations we show that the matrix 



where mi^ for A; = 5,6,7,8 are determinants of principal kxk 
submatrices of the matrix (IV.4I). It turns out that the fac- 



trB pa(0 «) 1 - Pa(0 



equals to 



/ ni 











-Pl5 











-Pl9\ 
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-P37 
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-Pl5 











^•55 











-P59 

















r66 





-P68 











-P73 











'"77 























-P86 





^"88 





\-P91 











-P95 











'"99 / 



where 



(V.4) 
P11+P22 + P33-PM, ^==1,2,3 

P44 + P55 + P66-PH, fc = 4,5,6 (V.5) 

P77 + P88 + P99-PM, fc = 7,8,9 



rkk- 



We compute that leading principal minors of the matrix (IV. 4b 
which can change the sign during the evolution, and find the 
following expressions 

ms =r22'"33'-44 ('"11 ^"55 - |pl5|^) 

me ='"22'"33'-44'"66 {rufSS - IPlsl^) 

mi =r22'"44'"66('"ll'"55- |Pl5p) ('"33'"77 - |P37p) 

mg ='-22'"44(ni'-55- |Pl5|^)('"33'-77-|p37p)('"66'"88- IPesI^) 

(V.6) 



tor rnrss — |pi5| is always positive, but as follows from the 
numerical analysis, r33r77 — |p37p as well as rggrgg — IpgsP 
change the sign during the evolution (see FIG. 2). Let to be 
the time at which the factor r66'"88 — |P68p changes the sign. 



We see that to > 0, so only after that time the matrix ( IV. 4b be- 
comes non - positive. It means that for t > to, the states pa(f ) 
are necessarily distillable. One can check that to > tN, and we 
see that the initial bound entangled state ( IIV.3I ) evolves in the 
remarkable way: for all f < f^r it is PPT, for t^ < t < to il is 
NPPT but a priori can be non-distillable and only after to it 
becomes distillable. To show this in the explicit way, let us 
introduce the measure of the violation of reduction criterion, 
defined as 

A^,ed(p)=max(0,-?ij^f„) (V.7) 

where XJ^ is the minimal eigenvalues of the matrix 

Pied = trg p (g) 1 - p 

The quantity ( IV.7I ) can be called the reduction negativity of 
the state p. For the bound entangled initial state ( IIV.3I ) the 
evolution of negativity and reduction negativity is given be- 
low (FIG. 3). So we observe in the system the phenomenon of 
delayed sudden birth of distillable entanglement. The phys- 
ical reason for the appearance of this phenomenon can be 
explained as follows. During the time evolution of the sys- 
tem, the process of the photon exchange produces coherence 
between the states |1a) <E) \3b) and |3a) <E) |1b), so the value 
of IP37I starts to grow. Similarly, the same process causes 
the production of coherence between states \2a) <E) |3b) and 
|3a) (8) |2b), so Ipesl also grows. Notice that non - zero value of 
|p37| or |p6g| is necessary for the possibility of creation of dis- 
tillable entanglement. But this condition is not sufficient since 
the populations of states of two-atomic system also evolve in 
time. We see from the formula ( |V.6b that there is a threshold 
for the reduction negativity (IV.7I) at which the system becomes 
distillable. 



0.020 




FIG. 3: The time evolution negativity A'o (dotted line) and reduction 
negativity A'jed (solid line) for the initial state ( IIV.3t with a = 3.6 and 
R = 0.2X. 



The numerical value of to depends on the choice of the pa- 
rameter a and the interatomic distance R. For the initial state 
with a = 3.6 and the distance R — 0.2 A-, we obtain foyw 0.78 
whereas InJ^ 0.49. 



the common vacuum. In the case of small (compared to the 
radiation wavelength) separation between the atoms, the sys- 
tem has nontrivial asymptotic states which can be entangled 
even if the initial states were separable. For the large class of 
separable initial states the asymptotic states are not only en- 
tangled but also distillable. The same is true for some class of 
bound entangled initial states. Thus we have shown that the 
dynamics of the system can transform bound entanglement 
into the free distillable entanglement of stationary states. For 
the atoms separated by larger distances only some transient 
entanglement can exist but still the dynamical generation of 
entanglement is possible. We have shown that this happens 
also for the class of bound entangled initial states. Moreover 
we have demonstrated that such states evolve in a very pecu- 
liar way: they almost immediately disentangle after the atoms 
begin to interact with the vacuum, then for some finite period 
of time there is no entanglement and suddenly at some time 
the entanglement starts to build up. But this entanglement a 
priori can be nondistillable. We have analysed this problem 
using the reduction criterion of separability and found that the 
free entanglement surely appears in the system after some ad- 
ditional period of time. 



VI. CONCLUSIONS 

We have studied the dynamics of entanglement in the sys- 
tem of three-level atoms in the V configuration, coupled to 
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